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superconductive state. These results may be compared with the experimentally de- 
duced self-energy from ARPES experiments. 
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1. Introduction 

One of the leading contenders for the <i-wave pairing mechanism of cuprate 
superconductors is the spin fluctuations. Superconductivity (SC) mediated by the 
spin fluctuations has a long history. [U EJ |3] In support of the spin fluctuation 
mechanism, Scalapino,[T] noticing the commonalities among the heavy fermion, cuprate, 
and Fe superconductors, argued that (a) Their chemical and structural makeup, 
their phase diagrams, and the observation of a neutron scattering spin resonance 
in the superconducting phase support the notion that they form a related class of 
superconducting materials, (b) A number of their observed properties are described 
by Hubbard-like models, (c) Numerical studies of the effective pairing interaction in 
the Hubbard-like models find unconventional pairing mediated by an S = 1 particle- 
hole channel. He proposed that spin fluctuation mediated pairing provides the common 
thread which is responsible for superconductivity in all of these materials. 

Along the same line have many works been published. Noteworthy is the work by 
Damn et a/. [I] They measured the spin susceptibility from inelastic neutron scattering 
(INS) experiments on YBa2Cu306.6 and used it as the effective interaction (the 
Eliashberg function) between electrons to calculate the diagonal self-energy from the 
Eliashberg equation. The calculated spectral function produced similar results as the 
measured angle resolved photoemission spectroscopy (ARPES) intensity from the same 
YBa 2 Cu306.6 crystal. They claimed that a self-consistent description of ARPES and 
INS can be obtained within the Eliashberg formalism for the cuprates (by adjusting a 
single parameter, the fermion-spin coupling strength.) Their work, however, did not 
take full consideration of the frequency and momentum dependence of the diagonal 
and off-diagonal self-energy. This point is crucial in that, for instance, the momentum 
dependence of the peak position of the self-energy which is one of the ongoing discussions 
in the field can only be addressed by calculations without assuming ad hoc momentum 
dependence. See the remarks in Sec. 5 below. 

Here, we revisit this spin fluctuation scenario by computing the angle, i.e., the 
momentum direction in the Brillouin zone (BZ), and frequency dependence of the 
diagonal, E(k, cj), and off-diagonal self-energy, 0(k, us). The cause of the angle and 
frequency dependence can provide an important clue about the pairing mechanism. [51 E] 
The diagonal self-energy is also called normal self-energy ("normal" here means the 
particle-hole channel and should not be confused with the "normal" as in the normal 
state meaning above T c ), and off-diagonal self-energy is also called anomalous or pairing 
self-energy. We employ the phenomenological fermion-spin coupling. (2j H] 



where a is the coupling strength of the dimension of energy, and c and a are the fermion 
and spin operators. Although the Eliashberg formalism is not firmly established for spin 
fluctuation mediated superconductivity, perhaps, our resort is that the ratio Xu c /Ef is 
< 1. A is the dimensionless coupling constant, u c the cutoff of the spin fluctuation 
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frequency, and Ep is the Fermi energy. Also, Millis argued in Ref. [TJ, in justifying 
the numerical Eliashberg approaches, that the (i-wave superconductivity induced by 
antiferromagnetic (AF) spin fluctuations is due essentially to high-energy part and not 
to the strong low-lying AF fluctuations producing the mass enhancement and scattering. 

For the Eliashberg function, a 2 F(q, u), we take like Dahm et al. the imaginary part 
of the spin susceptibility, Xspfaw), measured from INS. High quality INS data require 
large size single crystals and the INS data over wide momentum and energy range are 
mainly from YBCO or LSCO compounds. A functional form of the spin susceptibility 
obtained by fitting the INS data is given in the literature for optimally doped (OP) 
La2- 3; Sr 3 ;Cu04 (x = 0.16, T c = 38.5 K) by Vignolle et aZ.[H] and for overdoped (OV), 
(x = 0.22, T c = 26 K),[9] and underdoped (UD) Las^Sr^CuC^ (x = 0.08, T c = 22 
K)[10j by Lipscombe et a/.[llj The diagonal and off-diagonal self-energy and the quasi- 
particle (qp) energy shift, X(k,u), are computed self-consistently for OP, OV, and UD 
La 2 _ a; Sr 3 .Cu04 using the INS measured spin susceptibility. See Eqs. (fl2"l) and (J4]) below. 

In the following section 2, we present the Eliashberg formalism used to calculate 
the self-energy from the given spin susceptibility spectrum, x sp (q, cj). Some preliminary 
analysis for the energy scales of the self-energy is given in section 3 before presenting 
results of numerical calculations. In section 4, the numerical results are presented for 
OP, OV, and UD La 2 _ a: Sr a .Cu04 focusing on the angle dependence of the position and 
intensity of the peaks in the self-energy. There are two sources for the peaks in the self- 
energy in SC state: the peaks in the density of states (DOS) and the spin susceptibility. 
We will discuss how the two between them show up in the self-energy. The summary 
and outlooks will follow in section 5. 



2. Formalism 



The ci-wave Eliashberg equation is given by[T2l fT3] 

/OO /"OO 
de / de'Sj2 A s(k',e)a 2 F( + \k,k',e'), 
OO J —OO ^7 

/oo poo 
de / de'S^A x (k , ,e)a 2 F«(k,k , ,e'), 
oo J — oo 

/oo poo 
de / de , ,S^A (k',e)a 2 F ( - ) (k,k , ,e'), 
OO J —OO j^/ 

S - (2, 

e + e — oo — id 

where f(uj) and n(u) are the Fermi and Bose distribution functions, respectively. E is 
the symmetric part of the diagonal self-energy S, X the shift of the qp dispersion, and 
is the off-diagonal self-energy. The diagonal and off-diagonal Eliashberg functions, 
a 2 F^ and a 2 F^~\ are given by 

a 2 F«(k, k', e') = o4(k, k')F ch (k - k\ e') + a%(k, k')F sp (k - k' , e'), 
« 2 F(")(k, k', e') = o4(k, k')F c ,(k - k', e') - a% (k, k')F sp {k - k', e'), (3) 
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Figure 1. The Fermi surface of OP LSCO (black line) and OV LSCO (red line) given 
by Eq. ((TTJ) . The OP and UD LSCO have an hole-like and OV LSCO has a electron-like 
Fermi surface, respectively, as can be seen from the figure. The tilt angle 9 is given 
with respect to the nodal cut centered at (tt,tt) for OP and UD, and at (0,0) for OV 
LSCO. 



where the subscripts ch and sp refer to the channels by which the bosonic modes 
transform. For instance, the charge fluctuations belong to the ch channel, and the 
spin and current fluctuations to the sp channel. The various spectral functions are 
given by 

A s (k,u) = Im 



A x (k, oj) = Im 



7T W 2 — Y 2 - 



2 



7T W 2 ~Y 

M^^-'-Irn ^l^ , (4) 

where £(k) is the bare dispersion, 

W(k, u)=u- S(k, u) = tuZ(k, w), 

r(k,u;) = £(k)+X(k,u;), (5) 

and Z(k,u) is the renormalization function that appears, for example, in the gap 
function 

A(k,u)=<p(k,u)/Z(k,u). (6) 
The 4x4 matrix self-energy may be written as [14] 

S(k, u) = E(k, u)t + X(k, u)r 3 + <p(k, a;)r 2 a 2 , (7) 
where and <7j are the Pauli matrices in the particle-hole and spin space, respectively. 
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The diagonal self-energy is given by 

E(k,u) = E(k,u)+X(k,u), (8) 
and the diagonal spectral function measured by the ARPES is 

A(k,u) = A s (k,u) + A x (k,u). (9) 

The symmetry of the self-energy is as follows. 

E(k,w) = -E*(k,-w), 
X(k,u) = X*(k,-u), 

0(k,w) =0*(k,-w). (10) 

We choose the following tight-binding band as the bare dispersion for the doped 
La 2 _ x Sr :r Cu04. [15] 



£ (k) = — 2t(cos k x a + cos k y a) + At' cos k x a cos fc y a 

— 2t"(cos 2k x a + cos2k y a) — /x. (11) 

The tight-binding parameters are t = 0.25 eV, £' = 0.15 t, t" = 0.5 t', and the chemical 
potential /i = -0.19 eV for OP, and t = 0.25 eV, t' = 0.13 t, t" = 0.5 f , /i = -0.22 eV for 
OV, and t = 0.25 eV, f = 0.17 t, t" = 0.5 t', /i = -0.15 eV for UD Las^Sr^CuC^. The 
UD and OP La2_ x Sr x Cu04 have a hole-like and OV has an electron- like FS as shown 
in Fig. 1. To denote the momentum direction in the BZ in two dimensions we use the 
tilt angle 6 with respect to the nodal cut, that is, the diagonal cut along (0,0) — (tt,tt) 
line, centered at (71*, 7r) for the hole like FS and at (0,0) for electron like FS as indicated 
in Fig. 1. 

The spin fluctuation mechanism in this formulation means that we take F c h(q, u>) = 
and F sp (q,oj) as the spin susceptibility Xsp(Q) w ) measured by INS in Eq. (j3J). Then 
the imaginary parts of the self-energy may be rewritten from Eq. fl2]) as 

/oo 
de' [f(u - e') + n(-e')] ^ A s (k', u> - e')a 2 Xsp (k - k', 
-00 k , 

/oo 
de' [f(cu - e') + n(-e')} ^ A x (k', u - e')a 2 Xsp (k - k' 
OO y 

/OO 
de' [f(u - e') + n{-e')\ £ ^(k', u - e')« 2 Xs P (k - k', 
00 k , 

The real parts were calculated from the imaginary parts using the Kramers-Kronig (KK) 
relation. Eqs. ffl2|) and (j3J) were solved self-consistently via iterations. 

The coupling strength a was chosen such that it reproduces the experimentally 
measured gap amplitude A of L^-xSr^CuO^pi)] The coupling strength will be given 
in terms of the dimensionless coupling constant A below, 

A = f duN{u f Xs ^\ (13) 
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where N(u) is the density of states (DOS) and Xs P (u) is the local spin susceptibility 
given by 

= JW . (14) 

J d( l 

It is a measure of the density of spin excitation for a given energy. The gap function is 
determined by 

" 0(k,cu) " 
Z(k,u)_ 

and A was determined from the DOS peak position. 

The measured spin susceptibility x sp (q, u;) from INS was fitted by the form 



A(k) = Re 



(15) 

w=A(k) 



K 4 M 



•2| 



u) + fl(q)] 



2 • 



[{h-l)* + {k-iy-P] 2 + \{h-i)*(k-£ 



2 



i2(q) = " " ' V " " ^ ' " " - , (16) 

where the inplane wave- vector is written in the reciprocal lattice as q = ha* + kb* +/c*, 
k(u) = l/£ the inverse correlation length, the incommensurability S(u) specifies the 
position of the four peaks, and X(u) controls the shape of the pattern. A = 4 corresponds 
to four distinct peaks and A = corresponds to a pattern with circular symmetry. These 
fitting parameters were given in the references [HI EE HDl fTTj . Notice, however, that in 
these references the authors used the wave-vectors in the reciprocal lattice unit such 
that their 1/2, for example, corresponds to 7r of this paper. 

The measured local susceptibility may be decomposed into three parts; a low 
frequency incommensurate (IC) peak centered around = (7r±5, 7r) and the symmetry 
related points, a commensurate (CM) peak at Q = (±7r, ±7t), and a broad high 
frequency feature. The IC peak is around 18, 15, and 15 meV for OP, OV, and UD 
samples, respectively. The CM peak is around 50 meV for OP and 45 meV for UD 
samples, but is missing for OV La 2 _ x Sr x Cu04. On the other hand, the high frequency 
feature persisting up to measurable energy is common for all samples. The cutoff energy 
bj c of the susceptibility spectrum was taken to be 0.3 eV. This is the upper limit of the 
spin wave spectrum|17j of around 2J. 

The k' summation in Eq. ( fl2l) was performed by using the 2D fast Fourier transform 
(FFT) between the momentum and real space using the convolution relation 

e* k ' r E F ( k ' " k ) G ( k ') = F ^ G ^ ( 17 ) 

k k' 

on a 2 8 x 2 8 mesh of the first quadrant of BZ. No assumption about the k and u 
dependence nor a separable form of the diagonal and off-diagonal self-energy was made 
in the calculations. Self-consistency is reached in a couple of tens of iterations. 



3. Preliminary Analysis 



Before presenting our results of the angle and frequency dependence of the self-energy, 
it will be useful to consider some simple cases. Let us first consider the Einstein model 
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of frequency u b of the coupled boson. 

« 2 F(k, k', e') = a 2 (k, k') [5{e' - u b ) - 6{e' + u b )] . (18) 
Then the imaginary part of the diagonal self-energy from Eq. ( fl2|) is 

S 2 (k, lu) = ix {[f(u - u b ) + n(-uj b )} D(k, uj - uj b ) 

-[f(u + u b )+n(u b )]D(k,uj + uj b )}, (19) 

where 

D(k, uj) = J2 " 2 ( k ' k ')^(k', to). (20) 
k' 

In the low temperature limit of T — > 0, it is reduced to 

- S 2 (k, to) = ir [e(w - Lu b )D(k, uj-uj b ) + Q(-u - to b )D(k, u + u b )} , (21) 

where is the step function. The peaks of — S 2 (k, uS) for the negative (positive) u region 
are determined by those of D(k, u ± u>b). Depending on the range of k' summation of 
Eq. (1201 determined by a 2 (k, k'), either one peak (for k — > oo or 0) or two peaks (for 
intermediate h) may show up as discussed below. 

Consider two limits of this expression: First, for momentum independent coupling 
of a 2 (k, k') = a 2 . Then, we have 

{ixa 2 N(uo — Ub), for u > u b , 
ix a 2 N {us + u b ), for u < -u b , (22) 
0, otherwise. 

where 

N(u) = J2MV,uj) (23) 

k' 

is DOS. This clearly shows that the peaks of DOS at w = ±Ao in the SC state are 
shifted to ±(A + Ub) in — £ 2 (u;) because of the coupling to the boson of frequency of a;;, 
and that they are momentum independent. This case is relevant where the correlation 
length £ of a susceptibility peak is small (or, the inverse correlation length k — > oo) like 
the 50 meV CM peak of OP sample. See the angle independent peaks near \oj\ m 65 
meV in — E 2 and </>2 shown in Figs. 2(a) and 3(a). 

Another limit is where the coupling is a delta function like a 2 (k, k') = a 2 5(k.' — 
k — q), corresponding to k — > 0. Then, instead of Eq. ( |22l) . we have 

{7ra 2 A(k + q, to — u b ), for u > oj b , 
7ra 2 A(k + q, oj + u b ), for u < —u b , (24) 
0, otherwise. 

Because the spectral function A(k, uj) has a peak around uj ~ ±E(k) in SC state, where 



E(k) = Ve(k) + A*(k), (25) 

and £(k) = (£(k) + X)/Z is the renormalized dispersion, the peak of — S 2 (k, uj) occurs 
at uj pa ±(E(k + q) + uj b ) which is clearly momentum and band structure dependent. 



Angle and frequency dependence of self-energy from spin fluctuations 



8 



For intermediate values of k, the self-energy of Eq. f !24p is summed around k + q 
over the width k. Then, both energy scales of (A + Ub) and (E(k + q) + ojb) may appear 
in — E 2 (k, oj) from a single peak of x(q, oj). A modification is that the energy (A + w&) 
of £ — > now becomes momentum dependent (A'(k) + ojb) because a non-zero £ implies 
a momentum selection in the k' summation. A'(k) is an angle dependent energy of 
order A . This seems to be the case for the IC peaks as will be discussed below. 

Some complications arise from the band structure and momentum dependent 
coupling. [12] An interesting case is where the momentum sum covers the saddle point 
which usually occurs at (0,7r). This introduces another energy scale in — £ 2 (k, oj) 
from the van Hove singularity (VHS). It occurs at oj w sgn(^vHs)(EvHS + <^b), where 



Evhs = y£ 2 (0,7r) + Aq and sgn(f) = ±1 represents the sign of /, in addition to the 
two energy scales of the peaks discussed above. The shape of £(k, oj) are also modified 
by the impurity scatterings. (18] The VHS peak may be substantially suppressed by the 
coupling to boson spectrum and impurity scatterings. One should perform the self- 
consistent calculations to see their effects without misleading conclusion. The off-plane 
elastic impurities may induce interesting features of the self-energy in the SC state. 

In cases where the VHS peak is strongly suppressed and/or the Evhs and A are 
not well separated, the VHS feature may not clearly show up. The parameters of current 
LSCO calculations seem to belong to this case and we do not discuss the VHS features 
in the self-energy below. Also recall that the discussion so far is restricted to a sharp 
boson frequency of a single energy. A finite width in energy as well as in momentum 
space smoothens peak features in the self-energy. 

4. Numerical Results 

We now turn to the self-consistent numerical calculations using the experimentally 
measured spin susceptibility as the Eliahberg function for the OP, OV, and UD 
La2_ x Sr 3 .Cu04 as explained above. The angle 9 in the BZ was chosen with respect 
to the nodal direction as shown in Fig. 1. We wish to discuss the position and intensity 
of peaks in the absolute value of the imaginary part, — £ 2 (k, oj), and the real part of the 
self-energy, £i(k, oj), for oj < 0. The oj > region can not be probed by the ARPES 
with which we wish to compare our numerical results. For the off-diagonal self-energy, 
0i(k, oj) and 2 (k, oj) are, respectively, even and odd functions of oj, and will be shown 
in the oj > region. 

4-1. OP La 2 ~ x Sr x Cu 4 

First, we consider the OP La 2 _ x Sr x Cu04 with the doping concentration x = 0.16 and 
the critical temperature T c = 38.5 K. The spin susceptibility spectrum reported by 
Vignolle et al. has three parts; the IC peak near 18 meV, CM peak near 50 meV, and 
broad high frequency feature extended to 0.3 eV. The coupling constant was chosen 
such that A = 1.74 in the calculations to obtain the gap amplitude A = 17 meV in the 
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Figure 2. The calculated imaginary and real parts of the diagonal self-energy at the 
Fermi momentum of OP LSCO along several cuts in BZ in the SC state. The Vignollc 
spectrum was used for the F sp in the Eliashberg equation. The nodal direction is 
deg and the anti-nodal direction is 45 deg. 



Fig. 2(a) and (b) are the imaginary and real parts of the diagonal self-energy along 
several cuts perpendicular to the FS in the BZ in SC state. Fig. 3(a) and (b) show the 
imaginary and real parts of the off-diagonal self-energy, respectively. The two peaks in 
both the real and imaginary parts of the diagonal and off-diagonal self-energies reflect 
the two peaks in the spin susceptibility with slight complication as explained below. 




Figure 3. The calculated imaginary and real parts of the off-diagoanl self-energy of 
OP LSCO along several cuts perpedicular to the FS in the SC state. The angles are 
the same as the Fig. 2. 



For the Vignolle spectrum from OP La2-a;Sr x Cu04, the IC peak has an intermediate 
correlation length of £ ~ 4.1a, the CM peak has a small £ rs 0.8a, and the broad high 
frequency feature has £ ~ 0.7a. [8] From the discussion in the previous section, we may 
expect two peaks at ui ~ — (A'(k) + tow) and u[ ~ — (E(k + Qs) + Vic) from the IC 
peak, and one peak at u 2 ~ — (A +ucm) from CM peak. The u[ seems to overlap with 
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uj 2 , and two peaks show up in — E 2 (k, uj). The uj\ is expected to be angle-dependent and 
uj 2 to be ~ —67 meV because A = 17 and ujqm ~ 50 meV. This is what we obtained 
in numerical calculations as shown in Fig. 2(a). The exactly same argument holds for 
the off-diagonal self-energy, 2 (k, cu) as shown in Fig. 3(a). 

In order to understand the peak energy of the real parts of the self-energy, recall 
that the real and imaginary parts are related by the KK relation. It means that the peak 
energy of Xi(w) is shifted from that of — S 2 (w) by the width of the peak, that is, the peak 
energy of Si is expected at uj ~ — (A'(k) + ujic — Tjc) and uj ~ — (A + u>cm — ^cm), 
where T is the width of the peak. This is indeed what we obtained from numerical 
calculations. See the plots of £i(k, uj) and <fii(k, uj) as shown in Figs. 2(b) and 3(b), 
respectively. 

Now, we turn to the intensity of the peaks of the self-energy. From Eqs. ( I2TI) and 
( |20l) we see that — £ 2 (k, uj) is given by the sum over k' of a 2 (k,k')A(k' ,u + It 
means that there is large contribution from the k' sum to — X 2 (k, uj) if both k' — k rj 
and uj ~ — (E(k + Qs) + ujic) are satisfied. This is better satisfied near the anti-nodal 
region and — S 2 (k, uj) increases as the tilt angle increases for small \uj\. For large \uj\, 
however, either of the two conditions become ill satisfied and — S 2 (k, uj) is roughly angle 
independent. This is indeed what Fig. 2(a) shows. The IC peak in all four plots in Fig. 
2 and Fig. 3 is not highest at 45 deg but ~ 40 deg because of the incommensurability 5. 
On the other hand, the CM peak near 65 meV is highest at 45 deg as expected because 
the broad CM peak connects the anti-nodal regions most effectively. 

The angle dependence of the off-diagonal self-energy 0(#, uj) along several cuts 
is roughly d-wave like as shown in Fig. 3. The imaginary part looks like the local 
spin susceptibility Xsp( w ) with the suppressed high frequency part. The suppression 
of 4>2{0,uj) above ~ 0.1 eV shows that the high frequency part of the susceptibility 
does not contribute much to pairing because its broad momentum dependence is not 
very effective for d-wave pairing. The real part <pi(8,uj) increases as uj increases from 
and exhibits two peaks induced by the two peaks in the spin susceptibility and then 
decreases and makes a zero crossing near o; 2 . There, 2 (6 , ,o;) has a peak because of the 
KK relation. 

4.2. OV Lo2_ x Sr x Cu0 4 

The calculations were done for OV La 2 _ a: Sr a ;Cu04 as well. The difference from the 
OP case is that (a) the spin susceptibility spectrum does not have the CM component, 
and (b) the bare dispersion has smaller next nearest neighbor hopping amplitude of 
t'/t = 0.13 and has an electron-like FS. The IC peak at 15 meV of the spin suscetibility 
has the correlation length of £ ~ 2.7a and the broad high energy feature has £ ~ 0.6a. [10] 
The coupling constant was chosen such that A = 1.98 to obtain the gap amplitude 
A = 11 meV at T = in calculations. 

Fig. 4(a) and (b) are the imaginary and real parts of the diagonal self-energy along 
several cuts perpendicular to the Fermi surface in the BZ. The results may be understood 
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Figure 4. The calculated imaginary and real parts of the diagonal self-energy of OV 
LSCO along several cuts in BZ in the SC state. 




Figure 5. The calculated imaginary and real parts of the off-diagoanl self-energy at 
the Fermi momentum of OV LSCO in the SC state. 



like the OP case. The single IC peak in the spin susceptibility has an intermediate 
correlation length and induces two peaks in — E 2 (k, u). Because the correlation length 
is rather small, the angle dependence of the peak at U\ ~ —25 and U2 ~ —75 meV 
are weak. The peak at co± is the shift of the DOS peak because A + uic ~ 26 meV. 
The peak at cj 2 is from the — (^(k + Q^) + ujc) as was discussed in the OP case. We 
note that the VHS peak can not appear in the negative energy because OV LSCO has 
electron- like Fermi surface. 

Fig. 5(a) and (b) show the imaginary and real parts of the off-diagonal self-energy 
along several cuts in BZ. The angle dependence is roughly <i-wave like. The imaginary 
part of the off-diagonal self-energy $2(9,00) has a peak near u± ~ 25 meV and looks like 
the local spin susceptibility with the suppressed high frequency part as the OP case. 
The real part (pi(8,u) begins to increase as u increases from 0, has a peak and then 
decreases and makes a zero-crossing near u>i where the 4> 2 (9,u) has the IC peak. 
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4.3. UD La 2 ^ x Sr x CuO A 

The calculations were done for UD La 2 _ a Sr ;E Cu04 as well. The spin susceptibility 
spectrum reported by Lipscombe et aL[TQj [TT] for 8% doping La2- a; Sr, r Cu04 was used 
in the calculations. The spectrum Xsp{w) has three parts as the OP case; the IC peak 
near 15 meV, CM peak near 45 meV, and a broad high frequency feature extrapolated 
upto 0.38 eV. The IC and CM peaks have the correlation lengths of £ ~ 2.45a and 0.8a, 
respectively. [11 J The coupling constant was chosen such that A = 1.58 to obtain the gap 
amplitude A = 17 meV at T = in the calculations. |16] 




Figure 6. (a), (b) The calculated imaginary and real parts of the diagonal self-energy 
of UD LSCO along several cuts in BZ in the SC state. 



Fig. 6(a) and (b) are the imaginary and real parts of the diagonal self-energy along 
several cuts perpendicular to the Fermi surface in the BZ. The results may be understood 
like the OP and OV cases. Because the peaks of the susceptibility are broader in 
frequency than OP and OV materials, the peaks in the self-energy are not as sharp as 
the OP and OV cases. 




Figure 7. (a), (b) The calculated imaginary and real parts of the off-diagoanl self- 
energy. 
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Fig. 7(a) and (b) show the imaginary and real parts of the off-diagonal self-energy 
along several cuts. The angle dependence is roughly d-wave like. The two peaks in 
the diagonal and off-diagonal self-energy are from the IC and CM peaks of the spin 
susceptibility. Their energy in the imaginary parts is expected at u)\ ~ A + ujc ~ 32 
and ~ A + ujcm ~ 62 meV. The peaks in the real parts are shifted by the width. 
This is what we obtained from the calculations as shown in Figs. 6 and 7. 

5. Summary and Concluding Remarks 

We have calculated the full momentum and frequency dependence of the self-energy 
by solving the Eliashberg equation using the measured spin susceptibility from inelastic 
neutron scattering experiments on optimally, overdoped, and underdoped La 2 _ a; Sr a .Cu04 
cuprates in the SC state. The real and imaginary parts of the diagonal and off-diagonal 
self-energy were presented for several cuts perpendicular to the Fermi surface for each 
doping concentration. 

The results of calculations were discussed in terms of the angle (i.e. the direction of 
momentum in the BZ) dependence of the peak position and intensity. First, the angle 
dependence of the peak intensity is that the spin fluctuation induced self-energy S(k, uS) 
is very anisotropic in the momentum space for small \u\. We can see from Fig. 2(a) 
that the absolute value of the imaginary part of the diagonal self-energy of OP LSCO 
increases by a factor of about 5 from nodal to anti-nodal directions below the CM peak 
of | a; | < 70 meV. The real part also changes roughly by the same factor as can be seen 
from the inset of Fig. 2(b). This is understandable because the spin susceptibility with 
the correlation length £ rs 4.1a must mean that the quasi-particle dynamics is very 
different for different directions in the BZ. Second, The angle dependence of the peak 
position is that the CM peak is angle independent at around \ufy\ ~ 70 meV and the IC 
peak position increases from « 25 to 37 meV as the angle changes from the nodal 
to antinodal direction for OP LSCO. 

As alluded in the introduction, this angle dependence of the peak energy may not 
be properly addressed in the approaches where a separable form of the off-diagonal self 
energy like 0(k, u) = cj)(oj)[cos(k x a) — cos(k y a)] or <f>(\c,u) = 4>{uj) cos(20),[TJJ [7J H] or a 
phenomenological form of the spin susceptibility is assumed. For example, the separable 
form of the off-diagonal self-energy may give misleading results with regard to the angle 
dependence of an energy scale because the angle dependence was built in by hands. No 
such assumptions were made in this work. 

The present results may be checked against the proposed spin fluctuation theory for 
the cuprate superconductivity. The most direct evidence of the spin fluctuation theory 
will be to detect the angle and frequency dependence of the diagonal and off-diagonal 
self-energy experimentally and compare it with the results presented here. This will 
be an extension of the McMillan- Rowell procedure of phonon superconductors [H] to 
(i-wave pairing. The experiment of choice for this purpose will be the ARPES because 
of its high momentum and frequency resolution capability. 
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Indeed, the self-energy has been deduced by performing the momentum distribution 
curve analysis of the ARPES intensity for Bi 2 Sr 2 CaCu20 8+( 5. For La 2 _ x Sr :E Cu04 crystals 
high quality ARPES data are not available for comparison though. The requirement of 
high resolution ARPES intensity data is being met only recently for B^S^CaC^Og+s 
using the Laser ARPES. [221 EH EE] 

If we compare the present calculations on LSCO with the ARPES analysis from 
BSCO, bearing the difference in mind, we notice that the peak intensity of the real 
part of the self-energy from BSCO does not change so much as the present calculations. 
From the nodal (9 = 0) to 30 deg, the peak height changes by less than a factor of 1.5. 
The angle dependence of the peak position is different too. From nodal to anti-nodal 
direction the peak position decreases in contrary to the present calculations based on 
the AF fluctuations. |13| l21~] Proper comparison, of course, must await high resolution 
ARPES data from L^-xSr^CuC^ materials which may be quantitatively checked against 
the present calculations. 
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